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Abstract 

First principles calculations based on density functional theory (DFT) are performed to study the 
electronic properties of Mg acceptors in A1N at the presence of oxygen. It is found that Mg and O 
tend to form complexes like Mg-O, Mg2-0, Mg3-0 and Mg4-0 which have activation energies about 
0.23 eV lower than that of Mg (except of the passive Mg-O). The lower activation energies originate 
from the extra states over valence band top of A1N induced by the passive Mg-O. By comparing to 
the well-established case of GaN, it is possible to fabricate Mg and O codoped A1N without MgO 
precipitate. These results suggest the possibility of achieving higher hole concentration in A1N by 
Mg and O codoping. 
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Ultraviolet light-emitting diodes (LEDs) with wavelengths less than 300 nm are of con- 
siderable technological significance. They are potential alternatives to the existing gas 
lasers and mercury lamps in applications such as disinfection, air and water purification 
and biomedicine where the latter encounter difficulties due to their high operating voltages, 
low efficiency, large size and toxicities. l|] They also promise high density optical data stor- 
age and high-resolution photolithography. [2|] Wide gap semiconductors, diamond (E g =5.5 
eV) and A1N (E g =6.2 eV) are the two most studied materials for such LEDs and A1N is 
particularly favored due to its direct gap band structure and subsequent high light-emission 
efficiency. Diamond as well as A1N based LEDs have been fabricated recently. 3], 4] How- 
ever, in view of their low efficiency and high operating voltages, further developments are 
still required to improve these two LEDs to the point where they can be used as devices, j^] 

Fabrication of a homo-structured LED requires both p- and n-type doping of a semicon- 
ductor. Unfortunately, the asymmetry dopability of wide gap semiconductors makes the 
fabrication of LEDs very difficult. The asymmetry dopability means that a wide gap semi- 
conductor can be either p- or ra-type doped, but not both. [6] A1N can be easily n-type doped 
with Si. However, p-type A1N is of great challenge. The most promising acceptor for A1N 
is Mg. Yet, the activation energy (Ea) of the Mg acceptor in A1N is 0.5 eV. ?], i| As the 
ratio of carrier concentration to impurity concentration follows exp(— EA/ksT), where ks 
is the Boltzmann constant and T is the temperature, only a very small fraction ~ 10~ 8 of 
Mg impurities is activated at room temperature. Since the upper limit of Mg concentra- 
tion in A1N is 2x~10 20 cm _3 ,{4], the corresponding hole concentration for Mg-doped A1N 
would not exceed ~10 12 cm -3 . This concentration, however, is still well below that for device 
applications (which requires a hole concentration at least 10 17 cm~ 3 ). 

Recently, codoping has been proposed and applied to overcome difficulties in p-type 
doping in some wide gap semiconductors such as GaN and ZnO.Q] Significant improvements 
in hole concentrations have been achieved in these two semiconductors. In this approach, 
p-type dopants (D) are incorporated into the semiconductor along with a small amount of 
reactive n-type impurities as codopants {CD). Then in the host semiconductors, complexes 
like D-CD,D 2 -CD, D 3 -CD and D^-CD will form and more often they have lower ionization 
energies than that of mono-dopant D (expect for D-CD which is passive). Thus the hole 
concentration may be greatly enhanced. By this approach, the hole concentration in ZnO was 
improved from ~10 13 cm -3 in N-mono-doped ZnO to ~10 17 cm~ 3 or ever higher in (Al,Ga or 
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In)-N codoped ZnO.[l0j, llll, ll2| Codoping also enhanced the hole concentration in Mg doped 



GaN. The conductivity of M g do ped GaN could be significantly enhanced by annealing in 
an environment with oxygen. 13J] This enhancement in conductivity was attributed to the 



decreased activation of Mg acceptors due to the inclusion of oxygen and a corres pon ding; 
order of magnitude increase of hole concentration (from ~10 17 cm~ 3 to ~10 18 cm -3 ). IjJ] 

It is therefore reasonable to believe that the hole concentration of Mg-doped A1N can be 
improved using a similar approach. In this letter, we investigate the effect of atomic oxygen 
on the activation energy of Mg acceptors in A1N by first principles electronic calculations 
based on density functional theory (DFT). We first study the possibility of the formation of 
Mg n -0 (here n ranges from 1 to 4; atomic O occupies N site and atomic Mg the neighboring 
Al site, respectively) complexes in A1N. All these complexes under study are assumed to be 
electrically neutral unless otherwise stated. This is reasonable since even at an activation 
energy as low as 0.20 eV, only a small portion (~10 -4 ) of the acceptors are in charged 
state. Then the activation energies of complex acceptor Mg n -0 is calculated. To explain 
the results obtained, the density of states (DOS) are given . Most of the calculations are 
done using a 3x3x2 supercell constructed from A1N unit cell using the plane- wave DFT 



code VASP. 
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16] The lattice constants of the supercell are kept fixed to avoid effects 



from spurious volume expansion. The T-centered 4x4x4 A;- mesh is used for irreducible 
Brillouin zone sampling;. The ion-electron interaction is described by Vanderbilt ultrasoft 
pseudopotentials[17|]with local density approximation (LDA) for the exchange-correlation 
potential. The electron wave function is expanded in plane waves with a cutoff energy of 
400 eV. These parameters ensure a convergence better than 1 meV for the total energy. In 
all the doped supercells, atomic coordinates are fully relaxed using the conjugate- gradient 



18] until the maximum force on a single atom is less than 0.03 eV/A. 



algorithm 

For A1N in wurtzite structure the calculated lattice constants are 3.08 A for a and 4.94 
A for c with the internal parameter being 0.382, in good agreement with experimental 
values. Based on this the 3x3x2 supercell (Al3 6 N 36 ) is constructed. Previous calculations 
on formation energies have shown that O occupies N site in A1N as a deep donor {19I I20I] and 
consequently in our supercell one N atom is replaced by an O atom. This O donor may act 
as an attraction center to single Mg atom and form complexes such as Mg-O, Mg 2 -0, Mg 3 -0 
and Mg4-0 with Mg atoms occupying the nearest Al sites to O and without destroying the 
lattice structure. To study whether single Mg atom will bind to Mg n -0 (n—0,1, 2 and 3, 
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respectively) complexes, we define the binding energy as the energy required to form the 
M gn+i-0 complex from well separated Mg dopant and Mg ra -0 complex: 

A (n) = E(Al 36 _ n ^Mg n+1 N 35 0) + E(Al 36 N 36 ) 

-E(Al 35 M gi N 36 ) - E(Al m - n Mg n N 35 0) (1) 

where E is the total energy of the system indicated in parentheses. A negative A^ n+1 ) 
suggests that the Mg n+1 -0 complex is energetically favorable and stable while a positive 
/\( n ) suggests that Mg n+ i-0 cannot form. The calculated A^ n ^ are summarized in Table [B 
As can be seen, single Mg atom will bind to single O atom for a large energy decrease of 5.304 
eV. This large energy decrease can be attributed to the passivation of the extra electron of 
O by Mg. The energy decrease from single Mg atom binding to Mg-0 complex is 0.626 eV, 
suggesting that Mg 2 -0 complex will form provided the amount of Mg exceeds that of O. If 
there are still extra Mg acceptors available then Mg 3 -0 will form with an energy decrease of 
0.412 eV . Formation of Mg4-0 from single Mg atom and Mg3-0 is less likely although also 
exothermic due to a much smaller energy decrease of 0.157 eV. The binding of single Mg 
atom to the O atom in Mg n -0 complexes can be attributed to the larger electronegativity 
of O than that of N. 

Now we calculate the activation energy of the Mg2-0 complex which is the most prob- 
able complex when the concentration ration between Mg and O is around 2. The method 



proposed by Van de Walle is applied for the activation calculation: [21] 



E a = E tot [D-] - E tot [D°] — E v — AV[D] + E corr (2) 

where E& is the activation energy of the defect (donor) D. E tot [D~] and E tot [D°] are the 
total energies of the supercell with defect in charged (-) and neutral (0) state, respectively. 
E v is the valence band maximum of the bulk semiconductor; AV[.D] is a correction term 
to align the reference potential in the charged defected supercell with that of the bulk and 
it is derived from the electrostatic potential difference between the bulk and that of the 
defected supercell far away from defect site. Using a 5x5x3 supercell a value of -0.15 eV 
is calculated for both AV[Mp] and AV[Mg 2 — O]. Calculation obtained using a larger 
6x6x4 supercell gives out does not result in any significant change in these values. E corr 
is a correction term for the use of T-included fc-mesh sampling for the hexagonal lattice. 
In practice this is derived from the energy difference between the highest occupied level at 



r-point and other special fc-points (averaged) in the supercell containing the neutral defect 
D. The calculated E corr is 0.26 eV for supercell containing only Mg defect and 0.03 eV for 
the supercell containing Mg 2 -0 complex. The calculated ionization of single Mg acceptor 
in A1N is 0.40 eV which is in agreement with previous DFT-LDA calculation (~ 0.45 eV 
22 1) and experiments (~ 0.5 eVp, S])- For the Mg 2 -0 complex, the calculated E a is only 



0.17 eV, 0.23 eV lower than that of single Mg acceptor. If we apply this decrease to the 
experimental value, then following exp(— EA/ksT), the hole concentration can be increased 
by at least a factor of 10 3 . This is a significant improvement, although the total carrier 
concentration is still below the value desired. Other clusters Mg3-0 and Mg4-0 have very 
close activation energies to that of Mg 2 -0. Thus incorporation of some amount of oxygen 
into Mg-doped A1N can improve the hole concentration. 

The decreased activation energy in the Mg acceptor after attaching to O-Mg complex 
can be understood from the density of states of the defected supercell as shown in Figure. HJ 
Although O-Mg complex is passive and cannot accept the host valence electrons, it induces 
extra fully occupied states right on the valence band maximum (VBM) as indicated by the 
DOS curve of A1N supercell containing Mg-0 complex. In the Mg-0 complex some electrons 
have higher energy than that of the host valence electrons. Thus Mg atoms binding to this 
complex may be activated by electrons from these complex states rather than from the host 
states. So the activation energy is decreased. 

However, for such a codoping the growth condition should be chosen carefully. This is 
because under thermal equilibrium the formation energies of Mg n O complexes will be rather 
high. The formation energy is given by: 

Ef[D{q)] = E tot [D(q)] - E tot [bulk] - 

+q(E f + E v ) (3) 

where the E^ ect (q) is the total energy of the supercell containing the defect, E tot {bulk) is 
the total energy of a similar supercell containing the pure crystal, rij is the number of atoms 
that is involved in the formation of the defect with /ij being the corresponding chemical 
potentials. Ef is the Fermi energy which is set to zero at the valence-band maximum 
E v . The chemical potentials depend on the experimental growth conditions, which can 
be either Al-rich or N-rich. Formation of A1N crystal under thermal equilibrium requires 
HAi+^N=^Ai[bulk}+lfi N [N2\+AH[AlN], where AH(AIN) is the formation enthalpy of A1N. 



In the N-rich condition which is preferred for incorporating Mg at Al sites, the upper limit 
of \x n is given by ^[N^], i.e., the energy of N in a N 2 molecular. AH[AIN] is calculated to 
be -3.58 eV. The formation energy of Mg on Al site at neutral state is W [Mgf^(0)]=1.76 
eV with the solubility limit imposed by 3(j J M g +2fJ>N = 3fiMglbulk]+2fi N lN 2 ]+ AH(Mg 3 N 2 ). 
However, if O2 is present in N 2 flux, to avoid the formation of MgO precipitate the upper 
limit of fiMg follows: 

^Mg + fJ'O = fJ>M 9 [bulk) + Ho[0 2 ] + AH (MgO) 

where AH(MgO)=-Q.69 eV from our calculation. Under this constrain the formation ener- 
gies of Mgn-0 complexes become rather high. The formation energy of Mg 2 -0 complex is 
calculated to be 4.81 eV, suggesting an ignorable concentration of Mg 2 -0 complex in A1N in 
a growth process which is close to thermal equilibrium. However, this solubility limit can be 
overcome by proper growth conditions, i.e. high growth temperature and high growth rate. 
For example, in the doping of GaN with Mg as acceptors by gas-source epitaxy method, the 
inclusion of oxygen did not result in the MgO precipitate in GaN while improved the hole 
concentration significantly, jjj], Q] By calculation, the formation energy of Mg 2 -0 complex 
in GaN is 6.31 eV, significantly higher than that of A1N. Thus in A1N, the solubility limit 
problem is less severe. This suggests the possibility of Mg and O codoping in A1N without 
MgO precipitates. 

To summarize, we have studied the electronic properties of Mg acceptors in A1N at the 
presence of oxygen by first principles calculations. Our calculations suggest the formation of 
of Mg ra -0 complexes and their lower activation energies compared to Mg. Compared to the 
well established case of GaN, the MgO precipitate problem can be overcome. Our results 
suggest that the hole concentration in AlN:Mg can be greatly enhanced by oxygen codoping. 
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TABLE I: Binding energy for single Mg acceptor to Mg n -0 complex . Unit: eV 



n 





1 


2 


3 


aw 


-5.304 


-0.626 


-0.412 


-0.157 
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FIG. 1: Total DOS of supercells with pure A1N [curve marked by A1N] and A1N containing one 
passive Mg-0 complex [curve marked by AlN(Mg-O)]. The bold arrow indicates the extra states on 
the top of the VBM of the pure A1N. The Fermi level of the pure A1N is indicated by the vertical 
dashed line. 
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